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We create Bose-Einstein condensates of Rb in a static 
magnetic trap with a superimposed blue-detuned ID optical 
lattice. By displacing the magnetic trap center we are able to 
control the condensate evolution. We observe a change in the 
frequency of the center-of-mass oscillation in the harmonic 
trapping potential, in analogy with an increase in effective 
mass. For fluid velocities greater than a local speed of sound, 
we observe the onset of dissipative processes up to full removal 
of the superfluid component. A parallel simulation study vi- 
sualizes the dynamics of the BEC and accounts for the main 
features of the observed behavior. 
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Bose-Einstein condensates (BEC) in dilute atomic 
gases are macroscopic quantum systems which can be 
manipulated by a variety of experimental techniques . 
The current development of such techniques is opening 
up a wealth of possibilities to explore new physics, e.g., 
in non-linear atom optics and to study various as- 
pects of superfluid behavior in the precisely controllable 
context of atomic physics ||. 

Atoms confined in a periodic potential share some 
properties with systems of electrons in crystals. Effects 
known from solid state physics, like Bloch oscillations and 
Wannier-Stark ladders, have been observed by exposing 
cold atoms to the dipole potential of far detuned opti- 
cal lattices Macroscopic quantum interference has 
been observed in an experiment on a BEC confined to 
the antinodes of a far detuned optical lattice || . Bragg 
diffraction from a condensate has been induced in mov- 
ing optical lattices Q. This has been used, e.g., as an 
atom-laser outcouplcr and as a tool for spectroscopy 
of the momentum in BEC's ||. Applications of BEC's in 
periodic potentials range from matter-wave transport || 
to interferometry j^] and quantum computing p0[ |. The 
question of the stability of the BEC during the evolution 
in optical potentials is crucial for these applications and 
has been addressed in theoretical works [O. 

In this Letter we report on some novel aspects of super- 
fluidity in BEC's by studying their center-of-mass oscil- 
lations inside the harmonic potential of a magnetic trap 
in presence of a one-dimensional (ID) optical lattice. We 
identify different dynamical regimes by varying the ini- 



tial displacement of the BEC from the bottom of the 
trap. For small displacements the BEC performs un- 
damped oscillations in the harmonic potential and feels 
the periodic potential only through a shift in the oscilla- 
tion frequency. At larger displacements we observe the 
onset of dissipative processes appearing through a damp- 
ing in the oscillations. We can describe the experimental 
results in terms of an inhomogencous superfluid having 
a density-dependent critical velocity. In parallel we re- 
port numerical studies of the Gross-Pitaevskii equation 
(GPE), which capture the main features of the observed 
dynamics. 

In our experimental setup |L2| we now produce BEC's 
of 87 Rb atoms in the (F=l,m^= —1) state. The funda- 
mental frequencies of our Ioffe-type magnetic trap are 
oj x = 2tt x 8.7Hz and u>±_ = 2ir x 90 Hz along the ax- 
ial and radial directions, respectively. The condensates 
are cigar-shaped with the long axis (the x-axis) oriented 
horizontally. With a number of atoms N — 4 x 10 5 , the 
typical dimensions (Thomas-Fermi radii) are R x = 55 fim 
and R± — 5.5 /im. 

We create a ID optical lattice by superimposing to the 
long axis of the magnetic trap a far detuned, retrore- 
flected laser beam with wavelength A. The waist of the 
beam is two orders of magnitude larger than the short 
condensate axis, and therefore the resulting dipole po- 
tential in the condensate region has the form V(r) = 
Vq cos 2 (27rx/A). With a blue detuning S = 2ir x 50 GHz 
from the Dl-linc at Ao ~ 795 nm and an intensity 
1=1 mW/mm 2 in the antinodes of the standing wave, 
the dipole potential height of an optical lattice well is 
Vo/kB — 270 nK fL3|j . The spontaneous scattering rate 
in the antinodes is T sp ~ 0.7 Hz at this detuning and 
intensity. 

To prepare the atomic cloud in the ground state of 
the combined magnetic trap and optical lattice we first 
perform evaporative cooling in the magnetic trap until 
we reach a temperature slightly above the critical tem- 
perature, T ~ 1.5 T c . Then we superimpose the optical 
lattice and continue with the evaporation process down 
to a temperature T ~ T c /2, where the thermal cloud is 
no longer observable. We have checked that the time at 
which we switch on the optical lattice does only affect 
the atom number according to the spontaneous scatter- 
ing rate, but does not influence the BEC dynamics as 
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long as the potential is switched on at temperature above 
T c . It is important to note that, for the dipole potential 
strengths and atom numbers which are used in the exper- 
iments, the modulated atomic density does not vanish in 
the antinodes of the lattice but reaches a minimum value 
which is significantly different from zero. 
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FIG. 1. Density distribution of a BEC in a harmonic trap 
with a superimposed optical lattice, from a numerical simula- 
tion of the 3D GPE for TV = 3 x 10 5 and V /k B = 270 nK. The 
inset shows an enlargement of the central region of the BEC. 
The envelope of the modulated density distribution follows 
the parabolic distribution in the harmonic trap. 

Figure [l] shows the density distribution along the x- 
axis, as obtained by numerical propagation of the 3D 
GPE in imaginary time [^4|. Here, the condensate spans 
about 250 lattice sites. In the experiment, the density 
modulation on the length scale of A/2 cannot be resolved, 
due to the limited resolution (~ 7 /im) of the absorption- 
imaging system. The modulation on a short length scale 
raises the chemical potential to the value /i/fes — 170 nK 
for N = 3 x 10 5 in the combined trap. Instead, in the 
purely magnetic trap we have fx/ks — 47 nK. 

In order to investigate the dynamics of the system 
we translate the magnetic trapping potential in the in- 
direction by a variable distance Ax ranging up to 300 ^m 
by changing the currents through the coils of the mag- 
netic trap. The translation takes a few milliseconds, 
which is short compared to the longitudinal oscillation 
period 2tt/u> x . Therefore, the BEC finds itself out of 
equilibrium and is subject to a potential gradient which 
forces it into motion. The presence of the magnetic trap- 
ping potential ensures that the atomic cloud maintains its 
high density (maximum density n max ~ 1.5 x 10 14 cm ). 
After an evolution time t ev in the displaced trap, both 
the magnetic trapping and the optical lattice are switched 
off simultaneously and the cloud is imaged after an ad- 
ditional free expansion of 26.5 ms. The imaging beam is 
horizontal and directed perpendicularly to the long con- 
densate axis. From the absorption image we deduce the 
center-of-mass motion and gain information on the dis- 
tortion of the BEC. 

In the absence of the optical lattice, the center-of-mass 
motion of the BEC in the displaced trap is an undamped 



oscillation with frequency w x = 2ir x 8.7 Hz and ampli- 
tude Ax, to which in the following we refer to as the 
"free oscillation" . After switching on the optical lattice 
we observe dynamics in different regimes. 




FIG. 2. Superfluid oscillations of a BEC in the presence 
of an optical lattice potential of height Vb/fcs ~ 270 nK 
(squares) and in a purely magnetic trap (triangles) , for initial 
displacement Aa; = (31 ± 3) /im. The lines give results from 
a numerical simulation of the ID GPE at the experimental 
parameters. 

For small displacements, Ax <« 50 /zm, the dynamics 
of the BEC resembles the free oscillation at the same am- 
plitude but with a significant shift in frequency. Figure || 
shows a comparison of free oscillations and oscillations 
with superimposed lattice for Ax = (31 ± 3) /im. For 
the lattice potential Vb/^s — 270 nK we find a shifted 
frequency uj* = 2ir x (8.0 ± 0.1) Hz. As can be seen 
from Fig. 0, this frequency shift is also reproduced in 
numerical simulations of the ID GPE using an explicit 
time-marching method Jl4||l5| ■ 

The frequency shift can be explained in terms of a 
renormalization of the atomic mass in the band states 
originating from the periodic potential. From the data 
in Fig. H we obtain an effective mass m* /to = (uj x /uj*) 2 — 
1.18 ± 0.02. Different from earlier experiments on cold 
atoms in an optical lattice under constant accelera- 
tion 0, in the present small-amplitude regime under 
harmonic forces we are only exploring the states near 
the Brillouin zone center. The above value of to* refers, 
therefore, to states near the bottom of the energy band. 

The essentially undamped oscillations of the BEC on 
the time scale of the experiments in the present small- 
amplitude regime is a manifestation of superfluid be- 
havior. The coherent condensate is being accelerated 
through band states as if it were a quasi-particle [fl6|| . 
Also, for small displacements we observe only marginal 
heating effects, i.e., the small thermal cloud of atoms 
accompanying the BEC can be fully accounted for by 
spontaneous scattering. 

However, the BEC enters a regime of dissipative dy- 
namics when we further increase the initial displacement 
Aa; and hence the velocity of the condensate. In Fig. || 
we report the measured ratio between the first oscilla- 
tion peak amplitude and the free-oscillation amplitude 
as a function of the trap displacement, together with the 
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values from the numerical simulation. At a displacement 
Ax ~ 50 /im when the maximum velocity attained by 
the condensate is v ~ 3mm/s, this ratio suddenly devi- 
ates from unity, indicating the insurgence of dissipation 
in the condensate motion. As shown in the inset of Fig. 
||, the subsequent dynamics is a damped oscillation of 
the center-of-mass at a greatly reduced frequency. The 
damping increases by further increasing the initial dis- 
placement, as is seen in the main body of Fig. this 
behavior being also displayed by the simulation data in 
the same Figure. 
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FIG. 3. Ratio of the first-peak amplitude of the oscilla- 
tion of the ensemble to the free-oscillation amplitude, Ai/A, 
as a function of initial displacement Ax, for the potential 
V /k B ~ 270 nK and atom number N ~ 3 x 10 5 . Circles: 
experimental data, triangles: results from ID numerical sim- 
ulation. The inset shows a full oscillation for a displacement 
Ax = 60 /jm with (stars) and without (squares) optical lat- 
tice. Here, lines are fits to the data. 



the planes of the lattice. This results in a friction force 
which damps the motion. 

Let us therefore enquire about the relative number of 
atoms in the superfLuid component of a BEC in a state 
of motion at a given velocity v, N s (v)/N say, with N 
being the number of atoms inside the harmonic trap in 
the absence of the optical potential. In order to measure 
this function and to deduce a maximum critical velocity 
v m ax, we have varied the displacement Ace and recorded 
the atomic distributions after a fixed evolution time t ev = 
40 ms. For low velocities, up to about 2 mm/s, the sample 
follows the position of a freely oscillating BEC, the ratio 
N s /N being a constant approximately equal to 0.7. This 
reduction below unity is merely due to the loss of atoms 
by spontaneous scattering of photons from the optical 
lattice beams during the preparation of the BEC. 

Upon increasing the velocity of the BEC, we observe 
a retardation of a part of the cloud, leading to a well 
detectable separation from the superfluid component af- 
ter free evolution. For velocities v ~ 4 mm/s we observe 
that only the central part of the fluid is moving without 
retardation, leading to a drastically changed aspect ra- 
tio with respect to the "unperturbed" BEC. The spatial 
separation from the thermal component allows a clear 
demonstration of the superfluid properties of inhomoge- 
neous Bose-Einstein condensates and a precise measure- 
ment of the critical velocity. The data for N s (v)/N in 
Fig. ^ show a dramatic depletion of the number of atoms 
in the superfluid component as the velocity of the fluid 
increases. 



With the onset of dissipative processes the condensate 
shape in the experiment becomes distorted and a much 
broader distribution, compatible with a thermal compo- 
nent appears. A thermal component is not allowed to 
arise in the numerical simulation, which is still based on 
the GPE. However, the density distribution of the con- 
densate in the simulation becomes fragmented and its 
phase is completely randomized. That is, the condensate 
in this regime breaks up into subsystems residing in an 
essentially independent manner inside the various wells 
of the periodic potential. 

Superfluidity can be expected to disappear when the 
velocity of flow is sufficient for the spontaneous emission 
of elementary excitations, as is the case for a homoge- 
neous Bose gas In a trapped condensate moving at 
sufficiently high velocities, emission of phonons and other 
excitations is favoured and the gas becomes heated |||l8| . 
The essentially ID dynamics of the present sample im- 
plies an important role for longitudinal phonon excita- 
tions in these processes, with a spectrum of critical ve- 
locities because of the inhomogeneity. In a simplified 
picture, the optical lattice can be viewed as a medium 
with a microscopic roughness, which leads to a velocity- 
dependent local compression of the gas moving through 




5 10 
Velocity (rim's) 



15 



FIG. 4. The fraction of atoms remaining in the undistorted 
part of the BEC, N 3 /N, as a function of the velocity reached 
during the evolution in the periodic potential. The line shows 
a fit to the data assuming a 3D parabolic density distribution 
and a critical velocity proportional to ^n{r). 

To model the breakdown of superfluidity in the inho- 
mogeneous density distribution of the trapped BEC in 
Fig. |], we first discuss the position dependence of the 
longitudinal sound velocity in the sample. In an inhomo- 
geneous system, as already pointed out by Andrews et 
al. jn| in discussing sound propagation in a magneti- 
cally trapped BEC, one may define a local speed of sound 
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c s (r). This is given by c s (r) — yj (n(r)/m)(6fi/6n) pO| , 
where within the Bogolubov approximation the stiffness 
constant Sfx/Sn would be equal to the coupling strength 
g. We have evaluated the appropriate stiffness constant 
by giving a longitudinal stretching (or squeezing) to our 
model system by a relative amount e of order 1% at 
fixed transverse profile. We find Sfi/(fie) ~ 0.7. Tak- 
ing /x/fcs ~ 170 nK and Sn ~ —ne where n is the av- 
erage density (n ~ 0An max ~ 6 x 10 13 cm~ 3 ), we find 
Cs,max — 5.2mm/s for the maximum value of the local 
sound velocity at the peak of the BEC density. This is in 
excellent agreement with the data in Fig. |], showing that 
complete destruction of the superfluid component occurs 

&f Umax 

~ 5mm/s. 

Assuming, therefore, that the critical velocity v c (r) for 
local destruction of the superfluid component in the inho- 
mogeneous condensate coincides with the local speed of 
sound c s (r), we have v c (r) oc yj n(r). As observed, super- 
fluidity breaks down first in the low-density regions. The 
envelope function of the density distribution of the BEC 
is an inverted parabola in 3D (see Fig. [I]) and hence, 
by integration over the high-density region, we get an 
equation for the relative number of atoms in the super- 
fluid part of the BEC for a given velocity v, N s (v)/N = 
[5/2 x (1 - v*/vl ax f/i - 3/2 x (1 - v*/v* max )W] • 
This expression implies that about 90% of the atomic 
probability density is localized in a region which remains 
superfluid up to velocities v ~ w maa; /2. The line in Fig. ^ 
shows that the above expression for N s (v) /N gives a very 
good account of the data, the fitted value of the maxi- 
mum velocity being v max — (5.3 ± 0.5) mm/s. 

In further experiments we have also observed indica- 
tions that the dissipation onset occurs at higher velocities 
for decreasing Vq and that the BEC propagates without 
dissipation in a regime of very low atom number. We plan 
to investigate these behaviors in detail in future work. 

In conclusion, we have investigated the dynamics of 
BEC's in a periodically modulated potential, both ex- 
perimentally and in numerical simulations. By measur- 
ing the effect of the periodic potential on the sloshing- 
mode oscillation inside the harmonic trap we have de- 
termined an average effective mass of the atoms in the 
condensate. The combined use of a periodic optical po- 
tential with the harmonic confinement has allowed us to 
observe novel features of superfluidity in an inhomoge- 
neous atomic BEC and to demonstrate a new technique 
for measuring a local density-dependent critical velocity. 

The results of this work are of importance for future ex- 
periments using periodic potentials for the manipulation 
of Bose-Einstein condensates and for the understanding 
of dissipative processes in coherent matter waves. The 
precise control of the parameters promises to be a pow- 
erful tool for a quantitative exploration of novel regimes 
occuring at different atom numbers or tunnel barrier 
heights. 
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